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Abstract 

We determine the analytic expression of the damping rates for very soft mov- 
ing quarks in an expansion to second order in powers of their momentum in 
the context of QCD at high temperature. The calculation is performed using 
the hard-thermal-loop-summed perturbation scheme. We describe the range 
of validity of the expansion and make a comparison with other calculations, 
particularly those using a magnetic mass as a shield from infrared sensitivity. 
We discuss the possible occurrence of infrared divergences in our results and 
argue that they are due to magnetic sensitivity. 

pacs: ll.lO.Wx 12.38.-t 12.38.Bx 12.38.Mh 

keywords: hard thermal loops, soft quark damping, infrared sensitivity. 
ENSK-TP-18 



I. INTRODUCTION 

One runs into difficulties when one applies standard loop expansion to gauge theories at 
high temperature T: physical quantities like the dispersion laws become gauge dependent. 
Early work on the QED plasma using a hydrodynamic approach is [ij, followed by 0, 
taking account of one-loop quantum effects. The work discusses to one-loop order the 
QCD polarization tensor at high temperature and quark density and determines the gluon 



*abada@wissal . dz 



1 



dispersion laws to lowest order in the coupling g. It shows that these dispersion laws are 
gauge invariant but the one-loop-order gluon damping rates in the long-wavelength limit are 
not. It also shows that while chromoelectric Debye screening does occur to lowest order in 
the one-loop calculation, chromomagnetic screening does not, a gauge invariant statement. 
The non-screening of chromomagnetic fields at lowest order is also discussed in The 
massless- quark spectrum to lowest order in the coupling is described in ^ and the full 
quasiparticle spectra to lowest order at high T for the whole momentum range are given in 
0. The quasiparticle spectra are also described in for gluons and |jTO| for quarks. 

The problem of gauge dependence of the damping rates has been emphasized in works 
in which the gluon damping rates, particularly at zero momentum, have been calculated 
to one-loop order in various gauges and schemes and different results have been obtained 



TT|. It was realized that the problem is related to the fact that at high temperature. 



higher- loop diagrams can contribute to lower orders in powers of the coupling |]T2[ . In other 
words, the standard loop expansion is no more an expansion in powers of g^. In a series of 
papers, Braaten and Pisarski developed a systematic method for an effective perturbative 
expansion that sums the so-called hard thermal loops (HTL) into effective propagators and, 
equally important for gauge theories, effective vertices [1T3|-[15||. Using this method, the 
transverse-gluon damping rate 7t(0) at zero momentum was shown to be finite, positive and 
independent of the gauge [|16| . Later, a generating- functional formalism for high-T QCD in 
the HTL approximation was developed [l^ and a relation to the eikonal of a Chern-Simons 
gauge theory was found |T^. From there, a hydrodynamic approach showed that the HTL 
approximation is essentially 'classical' [p!9|] . 

Once developed, the important question to answer is whether the HTL-summed pertur- 
bation is reliable for calculations in QCD at high temperature. If so, it would constitute an 
adequate framework for describing the properties of the quark-gluon plasma. Of particular 
interest is the question of infrared sensitivity in massless gauge theories, worsened at finite 
temperature by the presence of the Bose- Einstein distribution which behaves like 1/k for very 
small gluon energies k: quantities tend to diverge like powers of the infrared cutoff rather 



than logarithmically as is the case at zero temperature pO[- This infrared problem is prior 



to the advent of the HTL scheme. It is for example shown in that infrared (and mass) 
singularities do occur but cancel out in first-order radiative corrections to the production of 
lepton pairs in thermal (massless) QCD. It is therefore most interesting to see if the HTL 
scheme is of any help in this regard. More precisely, does the HTL-summed perturbation 
constitute a workable framework in which infrared divergences are cured consistently order 
by order in the coupling? 

It turns out that the HTL summation dresses the massless quarks and gluons allowing 
them to acquire thermal masses of order gT, rrif and rrig respectively P,[7|-|l0[]. This means 
that to this lowest order gT in effective perturbation, the infrared region is 'safe'. But as 
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recalled, static chromomagnetic fields do not screen at this lowest order, they are believed 
to do so at the next-to-leading order g'^T, the so-called magnetic scale 0-1^ • Therefore, the 
test for the infrared safeness of the HTL-summed expansion starts really in next-to-leading 
order calculations. 

Much work has been carried regarding the use of the HTL summation up to next-to- 
leading order. For example, it is shown in that the HTL summation cures the infrared 
logarithmic divergence in the production rate of hard thermal photons in high-temperature 
QCD with massless quarks. Also, the damping of fast moving fermions has been thoroughly 
investigated in |gj2|j2gj2|j2|] for hot QCD and |26|,g^ for hot QED. In [gl for instance, 
it is first performed a bare (i.e., not HTL-summed) two-loop calculation and by taking the 
fermion slightly off-shell, there occurs to this order cancellation of infrared, both electric 
and magnetic, singularities. It is then shown that the dominant on-shell graphs are those 



dictated by the HTL summation. Also, it is shown in that it is in fact possible to 
find the leading contribution g'^Tlnl/g to the damping rates for energetic (hard) quarks 
and gluons in high-temperature QCD without having recourse to the full machinery of the 
HTL-summation program, a leading contribution already obtained for quarks in |12|, see 
also p3 |. 

It seems therefore that the infrared problem can somewhat be brought under control 
when it comes to describing fast-moving (hard) quasiparticles. But how about slow-moving 
(soft) quarks and gluons, particularly those on shell? A first indication of the sensitivity of 
the HTL perturbation to slow moving particles can be found in where the production of 
non-thermalized soft real photons in HTL-summed perturbation in high-T QCD is discussed. 
It is argued that this scheme fails to screen mass singularities in that it is not able to yield 
a finite contribution to leading order to the production rate, a physical quantity. However, 
the divergences involved in |28] are collinear in nature and come from dressed vertices. In 
this regard, an improved action has been proposed in which incorporates an asymptotic 
mass moo that removes singularities coming from light-like external momenta. 

As to the importance of the magnetic sector and the infrared sensitivity of the HTL 



scheme at next-to-leading order calculations, this is well demonstrated in the works ||3(]H32 



Indeed, calculates in HTL-summed perturbation the non-abelian Debye mass at next- 
to-leading order from the static limit of the polarization tensor. ||3T| determines the same 
physical quantity at the same order in the same scheme, but from the correlator of two 
Polyakov loops, a gauge invariant quantity. Paper [^] discusses the more general problem 
of next-to-leading order non-abelian Debye screening in one-loop HTL-summed perturba- 
tion. It argues that since the magnetic sector is nonperturbative in essence, the perturbative 
next-to-leading-order results may not be reliable. This is explicitly shown in, for example, 
the strong dependence of the analytic structure of the inverse of the static longitudinal prop- 
agator on the infrared behavior of the transverse gluons where the results differ significantly 
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depending on whether we regularize the infrared sector by introducing a magnetic mas^ or 
not. One other interesting point discussed is the important role of the magnetic mass in can- 
celling the gauge dependent terms when obtaining Debye screening from the Polyakov-loop 
correlator. Finally, comparison with lattice simulations indicates that the magnetic-mass 
enhanced results are more compatible with the lattice ones, hence the importance of the 
magnetic sector. 

However, Debye screening is static in nature. It is therefore interesting to examine 
dynamic on-shell quantities in order to understand better the infrared behavior of hot gauge 
theories in HTL-summed perturbation for soft moving quasiparticles. It turns out that the 
damping rates at lowest order for such (very) soft moving quasiparticles are quite suitable. 
Indeed, the HTL quasiparticle self-energies are real and so no damping is manifest at order 
gT; it starts at precisely the magnetic scale g^T. Thus, to exhibit damping, one needs to 
add to the inverse propagators the next-to-leading order contributions to the self-energies 
which are, as dictated by the HTL-summed expansion, one-loop corrections with soft loop 
momenta, hence all propagators and vertices have to be HTL dressed. 



The first such calculation is paper [16| which determined the damping rate 7f(0) for 



transverse on-shell gluons with zero momentum and found: 

7t(0) = 0.088iV,/T, (1) 

where Nc is the number of colors. The analytic calculation of 7t (p) to order {p/rrigf' where p is 
the momentum of the very soft gluon was carried in |^3[. The zeroth order (|l]) was recovered 
and it was indicated that the coefficient of the order [p/mg)^ may carry infrared divergences. 



The infrared sensitivity of the on-shell gluon damping rates has been emphasized in p^j35 
where the damping rate 7i(0) for longitudinal gluons with zero momentum was determined 
to lowest order g^T and found to be different from 7t(0) and infrared divergent. This 
is to be contrasted with the fact that at zero momentum, there must be no difference 
between longitudinal and transverse gluons This statement is emphasized in |26| where 



a Slavnov- Taylor identity for the gluon polarization tensor in Coulomb gauge is derived and 
when applied to the next-to- leading order gluon self-energy, the equality 7z(0) = 7t(0) is 
obtained]^. 

The next step must be the discussion of the damping rates for very soft moving quarks 
at lowest order g^T . This is because quarks are also important in the structure of hot QCD. 



^Introduced as an infrared regulator, a point we come to later in section three. 

^It is assumed in |36] that the spatial next-to-leading order HTL-summed gluon self-energy is 
isotropic at zero momentum. Our explicit and direct calculations do not recover this isotropy. 



This issue will be addressed in detail in ||37|. 
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Since they too acquire a thermal mass mj at order gT and their damping rates start at the 
magnetic scale g^T, it is all but legitimate to inquire about their infrared sensitivity. There 
are already the two works BH] and PP[ which determined independently the damping rates 



7±(0) for quarks with zero-momentum and found: 

7±(0) = aoC;/T, (2) 
where Cf = [N^ — l)/2Nc and ao is a finite constant depending on Nc and Nf, the number of 



flavors. For example, for Nc = 3 and Nj = 2, we have ag = 0.111 .. ., ||38|j39|l . Result (g) for 
quarks resembles result (|l]) for gluons. It is therefore interesting to determine the damping 
rates for moving quarks, but with very soft momenta. In this article, we attempt to obtain 
to lowest order g^T an analytic expression for the damping rates 7±(p), where p = pp is 
the momentum of the quark, using only as ingredients what the HTL-summed expansion 
dictates. 

Obtaining a compact expression for 7±(p) would be ideal but hardly feasible technically. 
Rather, we attempt to obtain an expression for 7±(p) in powers of p/nif up to second order. 
This expansion is carried early on in the calculation. We will specify in section three its 
range of validity and argue that in order to get an explicit expression for the damping rates, 
manipulating otherwise is practically intractable. In this work, we describe in detail how we 
obtain the analytic expressions for the first three coefficients ao, oi and 02 involved in the 



expansion, see (^4), and we defer the numerical evaluation of these to future work [gO|. This 
is because it necessitates the extraction of the potentially infrared- divergent pieces from the 
finite contributions, something somewhat complicated. An additional complication comes 
from the presence of two soft masses, rrig and rrif, and the discussion necessitates working 
each case apart. The numerical evaluation also necessitates the handling of potential di- 
vergences coming from soft light-like loop momenta. Experience with the transverse gluon 
damping rate 7t(p) indicates that these latter divergences may ultimately be brought 
under control, but the infrared ones would most likely persist. One interesting aspect to 
mention is that the order p/rrif in 7±(p) does not vanish contrary to the gluonic case. Pre- 



liminary results ||4^ tend to indicate that there are no infrared divergences in the coefficient^ 



ai but they tend to appear in 02. On the one hand, this puts 'some water' into our argu- 
ments regarding the acceptability of our early expansion in powers of the very soft external 
momentum. On the other, recalling that the infrared divergences we find in the gluonic 
sector start also at order [p/nigf' p^|35| , ^ , it would be interesting to try to understand 



why this is so. 

Finally, it is useful to recall that the intensive use of the damping rates as a mean to 



We already know from |3^,39| that the first coefficient ag is safe. 
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probe the properties of finite-T gauge theories is also due to the fact that in general, calcu- 
lations in the full HTL-summed perturbation are quite difficult beyond lowest order. The 
damping rates are the simplest such non-trivial quantities to handle. Indeed, though they 
come from one-loop graphs with dressed propagators and vertices, they are defined through 
the imaginary part of the effective self-energies, something that simplifies significantly the 
calculation. The only attempt to correct in HTL-summed perturbation the quasiparticle 
spectra to order g'^T we are aware of is that of M . 



This paper is organized as follows. In the next section, we recall the definition of the 
quark damping rates and write them in the context of the HTL-summed perturbation. Their 
determination amounts to that of the imaginary part of the next-to-leading order quark self- 
energy which we carry in detail in section four. Section three is devoted to discussing the 
expansion of the effective self-energy in powers of the external momentum and to the manner 
with which we regularize the infrared region. A comparison with other computational and 
regularization schemes is also carried, most particularly those shifting the pole of the effective 
gluon propagators with a magnetic mass. The final results are presented and discussed in 
the last section. 



II. QUARK DAMPING RATES IN HTL-SUMMED PERTURBATION 

We use the imaginary-time formalism in which the euclidean momentum of the quark is 
= (po, p) such that = Po+P^ with the fermionic Matsubara frequency po = {2n+l)iTT, 
n an integer. Real-time amplitudes are obtained via the analytic continuation po = —iu + O^ 
where u is the energy of the quark. A momentum is said to be soft if both u and p are of 
order gT; it is said to be hard if one is or both are of order T. The t/iree-momentum p of 
the on-shell quark is said to be very soft if p is much smaller than gT, say of the order g^T 
and smaller. We follow closely the notation of [|T^ and the HTL results we quote in this 
section can all be found there, see also p^JT5|] . 

The effective propagator for the quark can be written as: 

* (P) = - [7+pA+ (P) + 7-pA_ (P)] , (3) 

where are the euclidean Dirac matrices, 7±p = (7° ±i'j.p) /2 and A± = {Dq ^ Ds) '^ 
with: 

Do = Wo -Qo — ; 
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p \ p ^ 



(4) 
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^ In — The poles of 



X 



where the quark thermal mass is mj = ^JCf/8gT and Qo (x) 
A±{—iu!,p) determine the dispersion lawsQ ci;-|-(p) to lowest order in g. For soft quarks, one 
has: 



uj^{p) = nif 



1 ± 
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± . . . 



(5) 



At this lowest order gT, ^±{p) are real and the quarks are not damped. To get the damp- 
ing rates to their lowest order, one has to include in the dispersion relations the contribution 
from the effective quark self-energy *S(P). Therefore, the inverse of the quark propagator 
becomes: 



A-p^ (P) = *A-p^ (P) - *S(P) . 



(6) 



The effective quark self-energy has also the decomposition *S = 7° *P'o + il-V*Ds where 
*Dq and *Ds are the two functions to be determined in HTL-summed perturbation. The 
inverse of the quark propagator is then: 



7° Po+ */^o)+n-Pps+ 



(7) 



The damping rates for quarks are 7±(p) = — Imr2-|- {p) where VL± are the poles of 
Ap'(— 2f2,p). Since the self-energy *S is (7-times smaller than *A^^, we have to lowest 
order: 



7± {P) 



(8) 



w=Lu±(p)+iO+ 



where /± = -Do T Ds, *f± = *Do T *Ds and stands for d/duj. Using the expressions in 
(I), it is easy to expand the denominator in the above relation in powers of p/rrif. One then 
obtains: 



7± (P) = \ 



2 p 
1± - 



2 

3 nif 9 ^ 



+ 



Im*/± {-i^,p) 



(9) 



We see that determining 7-1- (p) to lowest order in g amounts to calculating the imaginary 
part of the next-to-leading order quark self-energy. 



(+) for real quarks and ( — ) for 'plasminos' 



only thermally excited quasiparticles. 
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The HTL-summed perturbation [jT3|-0 dictates that the next-to-leading order quark 
self-energy is given in imaginary-time formalism by: 



*S(P) = *Si(P)+ *S2(P), (10) 

where we have: 

*Si (P) = -g^CfTi^^i, [*T^ (P, -Q; -K) {Q) T'^ (-P, Q- K) *A^, {K)] , (11) 

and: 

*S2 (P) = -'-g'^C fTi,,,, [t^- (P, -P; ir, -i^-) *A^, (iT)] . (12) 

K is the soft gluon loop momentum, Q = P ~ K and Tr = T / 7 — with A;o = 2mtT, 

ko (2vr)'^ 

a bosonic Matsubara frequency. The subscript "soft" means that only soft values of K are 
allowed in the integrals; hard values have dressed the propagators and vertices. Note that 
since the loop momentum K is soft, both propagators and vertices involved in (|^) must be 
dressed. 

The effective gluonic propagator *A^u (K) is taken in the strict Coulomb gauge where it 
has a simplified structure. It is given by *Aoo (K) = *Ai (K), *Aoi (K) = and *Aij (K) = 
(^6ij — kikj^ *At (K) with *A; and *At having the following expressions: 

where SIli{K) = 3mJ(3,(!|i) and mt{K) = f [Q3V^) - QiC^) -§]. Q.(f' ) is a Leg- 
endre function of the second kind and the gluon thermal mass mg=^Nc + Nf/2 gT/3. The 
effective (dressed) vertices *r intervening in (^) are of the form: 

T = r + (5r , (14) 

where F is the bare (tree) vertex and 6r is the corresponding hard thermal loop. The two 
effective vertices that enter the calculation of the effective self-energy (0) are the effective 
quark-gluon vertex: 



.r.(p.g;«).,. + ,„5/;^_|^. (15) 



where the second term is the hard thermal loop, and the effective two-gluons-quark-antiquark 
vertex: 
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Note that the bare two-gluons-quark-antiquark vertex is zero so that the corresponding 
effective vertex is just the hard thermal loop. In both (|15D and (|T^) , S = {i,s) and Qg is 
the solid angle of s. 

The task is to attempt to get an expression for the imaginary part of the effective quark 
self-energy *S (P). The 'natural' sequence of steps to follow is first to perform the angular 
integrations in the dressed vertices (^) and ([TB|). Next is to do the Matsubara sum in (PUD. 
Only then the continuation to real quark energies po = + 0"*" can be taken and the 
on-shell condition enforced. Last is to find a way to perform the integration over the gluon 
loop three-momentum k. However, given the complicated expressions we are faced with, it 
is practically very difficult to follow this sequence of operations. What we do in this work is 
first expand the effective self-energy in powers of the quark momentum p/rrif. This allows for 
an easy angular integration over Qg. Only then do we perform the Matsubara sum, this by 
using the spectral representation of the different quantities involved. The angular integration 
over Qk is subsequently done and the remaining integrals can be calculated numerically. We 
discuss this procedure in more detail in the next section. 



III. REGULARIZATION AND EXPANSION IN EXTERNAL MOMENTUM 



The expansion in powers of the external momentum of the HTL-summed next-to- leading 
order self-energies can be questioned from the outset in view of the fact that infrared di- 



vergences do appear in next-to- leading order physical quantities like the damping rates [35 



Are these divergences genuine or merely artifacts due to the method used? First, recall that 
we are considering only very soft external momenta. Note also that the HTL framework 
itself allows for an expansion of quantities in powers of the soft external momenta. For 
example, the gluon and quark on-shell energies uj{p) are obtained in the literature in the 
form of a series in powers of soft p E3|,|T3[ . The same is true for the residue and cut functions 



intervening in the spectral decomposition of the effective propagators |^. It is therefore 
legitimate to expect the perturbation built on hard thermal loops (these being considered as 
a zeroth order approximation) to be analytic in very soft p, and hence admit an expansion 
in powers of such momenta. Also, such an expansion is not proper to this work, it has been 
previously used in the literature, for example in . 

There is of course a distinction between the analyticity in p and that in g. For example, 
the standard loop expansion of QCD is in powers of whereas in HTL-summed pertur- 
bation, the expansion is in powers of ^/g^. The same is true for other theories like the 
prototype X^p'^ theory and QED 0. This may introduce a non-analyticity with respect 
to g in some quantities, but does not necessarily change drastically the analytic behavior 
of these quantities with respect to very soft p. Take for example the estimation of the soft 
gluon damping rates made in 
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lt,i{p) ^^ln^t^i,«(p), (17) 

where Vtj(p) are the corresponding group velocities. These rates are clearly non-analytic in 
small g, but perfectly analytic in p: they even tend to zero0 as p ^0. 

Regarding the damping rates, our starting position is that the quark-gluon plasma is to 



be a stable phase of hadronic matter, at least for very soft excitations ||8|. QCD at high 
temperature in the (lowest-order) HTL approximation is 'finite'. At next-to- leading order, 
HTL-summed perturbation yields finite and positive damping rates for zero-momentum on- 
shell quarks and transverse gluons. The stability criterion ensures that we must expect the 
damping rates to remain finite and positive for non-zero very soft momenta. This translates 
into expecting the damping rates to admit a series expansion in powers of these very soft 
external momenta. This of course does not rule out a possible loss of analyticity for larger 
values of p, even just soft values. That would simply indicate new physics to explore. But 
because of the stability criterion, the analyticity must be preserved for very soft momenta. 
This is one important check to use in order to discuss the consistency and completeness of 
a given calculational scheme like the HTL-summed perturbation. 

Expecting an infrared problem, we introduce an infrared cut-off r] > such that /g^°° dk 
in (pUD and (|12D is replaced by J^°° dk. The cut-off r] is fixed for the rest of the calculation. 
It is physically useful to see it as representing the magnetic scale g^T. This means that k is 
never smaller than rj. In other words, we are summing contributions from all soft momenta^ 
k but not the very soft ones, i.e., those smaller than t]. We always regard the external 
momentum p as smaller than rj. We are therefore always working in the kinematic region 
< p < f] < k. This allows for the expansion in powers of p of all quantities that are 
functions of g = |p — k|. This is true in particular for 1/QS and the effective propagators 
*A (Q). The expansion of 1/ PS does not pose a problem in itself. 

It is useful to emphasize once more that our calculation sums the contributions from only 
the soft integration momenta rj < k: the very soft momenta < k < rj are systematically 
excluded and the hard region is cut by the spectral densities [^. If the integration is not 
sensitive to the very soft region (the magnetic sector), then the subsequent limit r/ — > in 
the final result should be smooth. If on the contrary there is sensitivity, it would mean that 
important contributions from this sector may be 'missed' by the HTL-summed perturbation. 
This is the essence of our point. Indeed, recall that regarding the self-energies, the HTL 
scheme discusses the two scales T (hard) and gT (soft), whereas with g and T, one has a 



hierarchy of scales g^T with n a nonnegative integer. As a matter of fact, argues that 



non-acceptable limit as we will discuss shortly. 
^Hard momenta are already summed in the hard thermal loops. 
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there may be scales between T and gT that play a significant role, and so n may not even 
be an integer. Since the HTL scheme, when built, does not consider effects like magnetic 
screening which (are believed to) arise nonperturbatively at order g^T and are not present 
at the HTL level, a perturbation built on the HTL summation may not be able to reproduce 
them. What it can do is to bring what contributes from the soft region to the very soft 
one. Therefore, excluding the very soft region from the fc-integration as we do may not be 
all unreasonable a thing to do. We think that the presence of magnetic sensitivity which 
manifests itself in infrared divergences indicates the very presence of these magnetic effects 
that the HTL scheme seems to be not able to accommodate. 

Let us compare our calculation of the damping rates with the estimation (0) mentioned 
above. This latter is different in many respects from the one we carry. It is obtained in the 
kinematic region where the loop momentum k is restricted to the very soft region (0 < < 
in our notation) whereas p is just soft, of order gT [0. In some sense, there, it is the very 
soft momenta that are integrated out; the soft ones are disregarded, something opposite to 
what we do. Result (|T^ cannot be carried to the very soft region p < in particular to 
the point p ^ for it will give zero (using the expressions of the group velocities at very 
soft momenta) whereas the damping rates there are finite. At the same time, our results 
can never be carried to the region p > t]. Clearly then, it should not be problematic if 
different analytic results are obtained. In fact, if really different results are obtained, which 
is the case, it only constitutes a further indication of the sensitivity of the HTL scheme to the 
magnetic sector. Equally interesting to note in the estimation (|T^) is that, in order to screen 
the divergent behavior at very soft momenta k, a regularization is used, which amounts to 
introducing a magnetic mass m^ag in the otherwise divergent propagators. It is clear that 
screening of chromomagnetic fields, if it occurs, is not necessarily going to manifest itself 
by a simple shift of the pole in the corresponding propagator by a momentum independent 



magnetic mass |32 



The presence of magnetic effects not handled by the HTL scheme is discussed in K 



It is argued there that for distances to order 1/T, ordinary perturbation (the standard 
loop expansion) is reliable. For distances to order 1/gT, the effective theory that screens 
static chromoelectric fields (the Braaten-Pisarski scheme) is reliable and it can be treated 
in perturbation. However, for distances to order 1/g'^T, one needs another effective theory 
which cannot be treated by perturbation (treated by lattice simulations for example). This 
last statement is emphasized by comparing the asymptotic behavior of the Polyakov-loop 
correlator determined from a magnetic lagrangian with that determined from an electric 
lagrangian with a put-by-hand magnetic mass mmag- With the magnetic lagrangian, an 
exponential decay governed by the lowest glueball state is obtained, an asymptotic behavior 
different from the one obtained from the magnetic-mass enhanced electric lagrangian. A 
comparison of these results with lattice simulations indicates that the glueball-state result 



11 



is more compatible with the lattice ones. One interesting inference one can draw from the 
above comments is that regularizing the infrared sector in HTL perturbation with a simple 
magnetic mass may not be the best description of magnetic effects, in particular if those are 
not incorporated in the scheme itself. 

It is important to stress that from a pure computational standpoint, matters are not 
straightforward if we defer in the self-energies the expansion in p after the angular integrals. 
Indeed, one has to deal with expressions quite complicated and involved, something of 

the sort T J2ko I J dils, J d^s2 ps,^f3s2QS2 *^ *^ ^^'>- Matsubara sum can be 
performed if the effective propagators, and are replaced by their respective spectral 
decomposition, see below. But this will bring in more than one energy denominator, which 
would compromise the straightforward extraction of the imaginary part of the effective self- 
energy. More serious a problem is the subsequent angular integration which will be very 
difficult, it not impossible, to perform [p^ . 

Finally, it turns out that for quarks, the expansion is in powers of p/rrif and not in 
(p/m/)^, as is the case for gluons (where m/ is replaced by m^), Preliminary results 

P0| tend to indicate that the second coefficient (that of p/ruf) is infrared safe together with 
the first one. This may suggest then that the expansion in powers of p is not sole to 'blame' 
for obtaining infrared divergent damping rates; other effects may be in play. 



IV. IMAGINARY PART OF ONE-LOOP HTL- SUMMED QUARK SELF-ENERGY 

Now we present the calculation of the HTL-summed next-to-leading order quark self- 
energy from which we extract the imaginary part. We first describe how we get an expression 
for Im *Si (P) defined in (|llD and then for Im *Zl2 (P) defined in (0). From now on, we 
take rrif = 1. This will simplify the final expressions we obtain. There remains another soft 

mass in the problem, m^, and so we define m = vrig/mf = K is easy to see 

that we always have m > 1. 



A. Calculation of Im (P) 

Using the structure of the fermion propagator (^) and that of the gluon propagator in 
the strict Coulomb gauge given just before (|T3p, we see that *Si (P) is composed of four 
terms: 

*Si(P) = ^ E T ^ y ^ [*r° (P, -Q; -K) A, (g) 7., T° (-P, Q; K) {K) 



e=± ko 



+ *r (P, -Q; -K) A, (Q) 7,, T^' (-P, Q- K) (5,, - * A^ {K) . (18) 
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The first two terms denoted *S£;(P), those with the longitudinal gluon propagator, are 
calculated separately from the two others denoted *T,^t{P)- We will illustrate the different 
steps of the calculation for *T,_i{P). Using the definition of the effective vertex (O) and 
making the change of integration variable K P — K, we have: 



— ry 

feo 



Sl 



dn 



7+fc - 



dQ, 2z^+707.k^+^7.k70 



4:71 



S2 



1 7-fc 



471 PS1KS1PS2KS2 



2PSKS 



A {K) *A, (Q) 



(19) 



The tilde over the sum sign indicates that fco is fermionic. 

Let us start with Ji, the term in (|19]) where there is one solid-angle integral over Vlg. This 
latter is carried in a reference frame where k is the principle axis (i.e., the '2;-axis'). The 
solid angle is then Vtg = {6,!f) such that k.s = cos^ and p.s = cos ■i/' cos ^ — sin'?/'sin6'sin(y9, 
where cos ip = k.p . Also, we have 7.S = 7'^ sin 6 cos + 7'^(sin ip cos 6 + cos ip sin 6 sin + 
j'^ {cos ip cos 6 — sin sin ^ sine/?), where {7'*} are the three spatially rotated Dirac matrices 
written in a reference frame where p is the principle axis and k in the {y, 2;)-plane. They 
are fixed in the integration over Qg. Performing all the (anti) commutations, we have: 



47r 



PSKS 



(20) 



In order to be able to perform with ease the above solid-angle integral, we use the expansion: 



1 
P5 



1 



p.s p.s 

1- — + 



Po 



(21) 



This expansion is valid in the region p < \ipo\, a condition always satisfied before analytic 
continuation and after. Before because po = (2n + 1) vrT and p ~ g^T. After because for 
very soft momenta, ipo = fnf + 0{p/mf) ~ gT, see (^. The solid-angle integral in ( pO] ) then 
reads: 



, r „7 (1 — cos6') -I- Z7.k — Z7.S 

dip / dO smO ; ; 

^Po -^0 Jo ikq + k cos 



X 



p p 
1 (cos-i/'cos^ — sin sin ^ sin ^■(cos-?/' cos^ — sin-i/'sin^siny^)^ -|- . . . 

Po 



The angular integrations are now straightforward and we obtain: 
d^k 1 



h 



7 



-1+ 1 + ^ o„.-^^^ 1+ , 

V k ipo \ k 



+ 



Po 



kl kl ( iko\ ^ \ 1 , 







x- - 




1 — -J^Qok 
k"^ 

_ _2. 

fc2 
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+ 1 + 



iko 
k 



+ 27. k 



Qok — 



px 



— Vofc H — 2 
/ Po 



k 



^iko iko 



+ 



Po 



X 



-Z7 



/3 



^3 
X ( 1 



+ 2 



A;2 



. ,z • , 1 iko px 

27 sm 1 — + — — 

k Ztpo 



^ko \ p 
k I tpo 



ikn 



kn 



X 



Ok 



+ 



Po 



'ikn 



X 



ikl 



3k k^ 
xA_ (K) *Ai (Q) , 



k^ 



Ok 



-x(l 



+ . . 



(22) 



where x = cosip and Qok stands for Qoiiko/k). 

The next step is to perform the integrals over the sohd angle of k in a reference frame 
where p is the principle axis. For this, it is most useful to develop all functions of g = |p — k| 
around k for (very) small p. The validity of these expansions is discussed in the previous 
section. In particular, here we need: 



(go, q) 



i-pxdu + ^^{^-^du + x'di\ + 



*Ai,t{qo,k) . 



(23) 



The solid angle of k is = (■?/;, a) and we have the relation: 



7 



7"*^ cos a — 7^ sin a; 



7' 



/2 



7"*^ sin a + 7^ cos a; 



7' 



/3 



7 



(24) 



where the {7'} are the (fixed) spatial Dirac matrices. Using (p3D, the integrations over if] 
and a become straightforward. We obtain: 



h 



A r+co h 



fco 



7 



1 + ^) e„. 



(go, k) 



(25) 



where A-k stands for A_{K). Note the introduction of the infrared cut-off rj. Note also 
that terms proportional to p do not vanish, contrary to what happens for gluons [p3| , p4 . 



The next step for Ji is to perform the Matsubara sum. This will be done after we get 
for the first term Jo in *S_;(P) (the one that involves no angular integrals) and the third 
term I2 (the one that involves two such integrals) expressions similar to (^). As for Jq, the 
calculation is simpler: only an integral over VLu using the expansion (|23D is needed. We get: 



14 



37r2T2 



ko 



dk k'^A_K 



7 



3 + P^T9k + 7;dl 



il^pdk + • • • 



'Aiiqo^k) . (26) 



As for /2, more work is needed. But we are fortunate here since Si is not 'coupled' to S2 so 
that each sohd-angle integral can be performed independently from the other. Each integral 
is thus performed along the lines shown for Ji, and so, there is no need to re-display the 
steps. After the two integrations are done, we multiply the two results, keeping terms to 
order only and taking care of the Dirac algebra. We obtain: 



7r2T2 



ko 



dk 



( ikn kn ' 

x*A/ {qo,k) , 



02 . ■ 3P 2 



1 



oP 



Wo 



iko 
k 



a- ] dk - a- { ^dk + -dl 



2-6- 



+ . . . 



ikn 



(27) 



where we have denoted for short Qe = 1 + e {l — Qok , £ = ±. We can now put together 
Jo, /i and I2 to get a first expression for *S„^(P). Since *S+i(P) is calculated in the same 
way and the only differences are mere signs, it is more economical to write the result for 
both terms in one single expression. We find: 



2 — r+oo r / 



2e 



L_ 2 

,2 u-}. 



ipok p^k"^ 



■ ■^P 

' 3 



edk + 



asi- h <9fc - 



e 



ipok \ipo 



1 

k 



1 

2' 



vdk + ^-dl ) + 



3 - 2 1 + 3£ 



ipnk \ p^ 
iko 



Pok^ 
e 1 



al[ — + dk 

Vpo 



iko 



-^ae + — ( 1 - e^as ] dk + eae [ T^k + 7:dl 



Q ^0 \ 2 



a. + 5 + 2e—^ - 3 



2e 



iko 



— 



\ 



k 



fc2 



a. - —a. 1 - £— a. d, 



Wo 



k 



+ ... 



'Ai (go, k) 



(28) 



.k ' 2' 

Now we are ready to perform the Matsubara sum over fermionic ko. We will need the 
spectral decomposition of A^{ko,k), *Ai{qo,k) and Qo{iko/k). They are worked out in 
,|T5[ and are given by: 

nl/T 



Ae{ko,k) 
At,i{ko, k) 
Qoiiko/k) 







l/T 
1/T 







dre'^''^ j dujp,{uj,k){l-n{uj))e- 
duj pt^i{uj, k) (1 + n{ijj)) e 
dTe'^°^ / duj po{uJ,k){l-n{uj))e~ 



1 



00 
+00 



00 
■ 00 



(29) 



n{uj) {n{ijj)) is the Bose-Einstein (Fermi-Dirac) distribution and the rho's are the spectral 
densities. Before replacing these above quantities, it is first necessary to rearrange terms in 
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(p8|) in a such a way that products of only two such functions appear. The reason behind is 
to ensure the appearance of only one energy denominator just before the extraction of the 
imaginary part, see below. For this purpose, we use the following easy-to-check relations: 



a^As = -ek [1 - {ik^ - ek) A^] ; 

a^Ae = —ek [a^ + ek {ik^ — ek) [1 — {ik^ — ek) A^]] . 



(30) 



After rearrangements, many terms happen to be real. Dropping these will yield: 

2 



Im *T.,i{P) 



Im TJ2 

ko 



dk k^ 



7 







1 + 



ik^ — ek 



A, + ^4 

vik 



-73 



^ 3 



ikn — ek ( ^ ikn — ek 
2——^ 1 H 



A,-£ 1 + 



ikQ — ek^ 



1 



2e 



Asdk - ^al I — + 9, 



Wo 



2p^k^ 



^iko — ek ^ {iko — eky 



2pl 



^ ikp - ek / ^ g^ifcp' 
Pok \ k ^ 



/ 1 2 ikn (iko — ek) ikn — ek 
\k p^Qk 



iko{ikQ — ek)^ ikQ — ek {iko — ek)'^' 



ipok 



ipok 



2ptk 



5 + 2e- 



ikn 



Pok 
iko 



Pok 



Aedk + 21^ + 



ipok 
iko — ek^ 



Wo 



k'^j'^ ^ \ ipok"^ Pok'^ , 



a'dk 



2plk 



'Aiiqo,k). (31) 



Here = — 1 = e (1 — Qok- Since iko appears in (|31|) only in the numerator of 
fractions, we can sum over it using the spectral decompositions (^9]) . At each time, we are 
left with two frequency integrals together with the one over k. Now we are allowed to take 
the real-energy analytic continuation ipo ^±{p) + ^0^. But just before, every has to 
be replaced with —1 except in the energy denominators which occur only once in each term, 
thanks to the rearrangements we made using (^Dj). The extraction of the imaginary part 
becomes straightforward if we use the relation 1/ (x + ^0+) = Pr (1/x) — iird^x). We obtain 
the following expression: 



Im *E,i{P) 



vrT 



dk 



+00 



duj 



duj' 



00 UJ 



-5 {uj±-u- u') [7° (A;2 {l + u- ekf p. 



{k - euj) po) + ^ [^27° (^ - ek) [k^ {1 + 00 - ek) p, + epo) 



+ iY [2k^ {00 -ek){l+uj- ek) p^ + ek^ {1 + u - ek) pA + {uj - ek) po {2e + dk) 
2 



3 



7 



3 2 1 

--e{uj-ek){k + 3euj) + {uj -ek)+-{uj- ekf (3k'^ + 2eujk + 3uj' 



+A; (3 + 4 (cj - ek) + k 



2 uj' 



2 

k'^ 

2uj {uj - ekfj pedk + — {I + uj - ekf p^dl 



- ek) (sP + 2eujk + 3u;^) po-^{uj- ek) {2uj - 1) podk + ^ {uj - ek) pod] 
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2A;2 {uj - ek) (2 + 3cj - 3eA;) p, + 2ek'^ {I + u - ek) {uj - ek) pA 
+ 2{uj-ek)po {3e + dk)]] + . . .] p'l . 



(32) 



Recall that we have set uif = 1. The notation is as follows: p^fl = Pe,o{^, k); p[ = pi{uj', k). 
In the above expression, we have used h{uj) ^ ^ and n{uj) ^ T/uj. This is because only soft 
values of u and u' are to contribute. The resulting integrals are to be performed numerically, 
but after the extraction of potential infrared divergences. 

It remains to calculate the two other contributions to Im* Si, those coming from trans- 
verse gluons in (p^ . The final result is the following: 



2 f+oo f+co f+oo fjrj' 

Im *S,i(P) = — / dk du —5{uj±-uj- u') 



7 



--{2k + er 



-^ {k'- -e{2k + e) {e - ) P. + ^ {eu^ + k) {k' - u') po 



+ 



P 



k'^ - +2(1+ ek) (k^ - cu^) + 26k + 1] p 



+ Y2^k + eu) {e - u') Po] + [(- {k' - ^')' - | {k' ^ ' 



[k^ + eujk 



+ 4cj(cj - ek) + 4£y - 2cj - 3 - pe + £ Q (/c^ - u?^ + 2ek (k 



u? 



k^ k^ 



+ 2ek + -)pA+ {2— + — {k + euj)) ik' - uj') po - -- {k + eoo)' {k - eu) pod, 



2P 



p2 
k 



1 



UJ 



UJ' 



lie 



UJ 



Uj'^ UJ^'' 



-e (u - ek)' i + s- l + e- +kiuj-ek) -k + —uj-- -3e- +- 



A;3 



^7e I3uj 475 



UJ" 



-iuj-ek) ------uj'-- -k - 



(Ae lOu 2e 



3k 3A;2 



A;3 



3k 3P 



2uj-^ 



A;3 



k k'^ k^ J 
A _2he u? 



+ 



lAe UJ 3 \ 
"saT ^ F ~ 2F/ 



-^(^^--^) {k^--'-^ek. 



k 



2e. 



ur Ae UJ 



+ - {u - ek) ( -k' - —kuj - 2uj' ] - -k + 2euj + tt + T ~ WT I Pe^' 



k 



2k 



+ ^ (- {uj' - +Aek {uj' - - Qk' - Aek + 2uj' - l) p,dl - ^ ((l + (l + 

{uj — ek) Po — e ( — 



UJ 



3 + 2.^-^ 



k P / 2F 



' ^ + 2.^;^-3-"^ 



A; A;2 



V2A;3 P 



tW[{-3 (k 



2 ^2 



2 4e 



(A;2 - cj2) (k'^ + £A;tu - cj^) +4 {to - eA;) (2cu + ek) 



- A (^ek + UJ -2£y^ - 5 - p^ + (^e (A;^ -u'^^ + 2{k + e) (k 
f 3 \ 1 



2 cu^ 



+ 



2A; + £ ) pe^fc 
(33) 
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This expression is quite long because these two terms are more involved. However, there are 
not new steps worth discussing in detail. 



B. Calculation of Im *T.2 (P) 



Now we turn to calculating the imaginary part of *S2 (P) which can be written from 



(|T^) and the structure of the gluon propagator in the strict Coulomb gauge as: 

d^k r dn, iS 



(P) 



ko 



{2^f J A-K KSPSQS 



A^i^)- 1-k.s^ *\{K) . (34) 



The steps to carry are similar to the ones we used for the previous contribution. But here 
we need the additional expansion: 



1 



QS iqQ — k.s 



p.s 



p.s 



iqo - k.s (^iq^ _ k.g 



+ . . . 



(35) 



The angular integrals over the solid angles fls and then Qk are done as usual. We obtain for 
the longitudinal gluon: 



'^21 (P) 



4 



7r2T2 



Po 



ko 

~~J?.Qok 



/ dk 



Po <J 



An' ^ 



2k 



Po 



+ 



iqok 



3ipo {ql + /c2) 3 {ql + k^) 



+ . 



^A^fccfc), (36) 



and for the transverse one: 



'^2t {P) 



IT' 



lj^2 



Po 



ko 
1 

,Po 



dk 



k 



l'-, {ql + k') Q'ok - ^I'X (- {ql + k') + ^qo] Q'ok 

Po ^ ' m \wo ^ ' j 

*^t{k^,k), (37) 



3.^0 + 3]^°^ + 



3po \Wo 
iqok 



3(g2 + P), 



where Qg^^stands for Qo {^^ ■ The sum over bosonic ko is now readily done if we add to the 
spectral representations (^) those of 1/ (gg + k'^) and 1/ {ql + /c^)^. We have: 



f-l/T . r+oo 

{ql + k"^) Jo J -00 



1 



(ql + k' 



l/T r+oo 

dre"^"^ 
J -00 



du e{iu)S{uj^ - (1 - h{uj)) e-' 
due{uo) 5^^\uo^ -k^){l-fi{uo))^ 



(38) 



with go fermionic. e(u;) is the sign function and S^^^uj"^ — k"^) stands for 5^^2(5 (u;^ — k"^). 
The extraction of the imaginary part is straightforward. We obtain for the longitudinal 
contribution: 
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Im *J:2i{P) = / dk duj 

TX 1 Jf) J —oo J —CO (jJ 



-^5{uO± — CO — to') 



±2l'kpo 



+ p 



7° (-f Po + - k') 



3 



+ . . . 



P'l, (39) 



and for the transverse one: 
4 r+oo 



A n+oo r+oo p+oo rj, / r 1 , 

Im *S2i(P) = / rfA; / du —5{u± - uj - uj') ^\ (k: 

TTI Jri J-oo J-oo UJ I k ^ 



3 L ' k 



+ p 

+ ... 

The final result we aim at is the sum of the six terms: 



Pt 



Im*S(P)= ^ Im*S,,(P)+ ^Im*S2i(P) 



(40) 



(41) 



where the different contributions are given in (|32|) , (133|) , (|39D and (^Ol) . To get the damping 
rates 7± (p), we use eq (P) where */± = *Dq =f *-Ds and *S = 7° *P)o + h-P *Ds. There are 
few more steps though. Indeed, note that the energy uj±{p) appearing in 5{uj± —u — uj') is a 
function of p, given in (|]) for small p. This means that for the terms in p^, the energy a;±(p) 
can be replaced by one (in units of m^) since we look for the damping rates up to order p^, 
but for the terms of order p, we have to expand 5{uj± — uj — u') to order p and for the terms 
of order zero to order p"^. A subsequent rearrangement is necessary. 



V. RESULTS AND CONCLUSION 



The damping rates are given in (Op. We find: 



7± (P) 



9'CfT 
Svr 



ao ± -ai + —a2 + ... 



where the coefficients Oj are given by the expressions: 
ao = 



00 n+oo n+oo 

dk I duj I — -fo {u, u'] k) 5 



V 



00 UJ 

fOO (^(^' 

— [ /i {uj, uj'; k) - /o (cu, J] k) dj\ S, 



roc /'+00 

ai = dk duj 

Jr] J— 00 J— 00 to' 

roo r+oo r+oo (j^f^j' ^ 

a2= dk du — f2{uj,u';k) -fi{uj,u';k)d^ +fo{uj,uj';k) 

J r] J —00 J —oo UJ 



(42) 



6, (43) 
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with S = S {1 — uj — uo'). The three functions fi (cj, uj'\ k) are given by the following expres- 
sions: 



2 2\ / 

^ ) PoPt 



(44) 



This expression is the one obtained in p^j39[] . /i and /2 are new. They read: 



/i {u, uj'; k) =Y, 2A;2 (-1 + k^ - 2ekuj + u'^) p^p'i + ( -'i + 2ek + 4k^ - k 



PePt 



- (2 + iek + 2A;2) u + i^— + 4 + 2ek + 2k'^j + 2uj^ - i^— + 1 
+ek\l- ek +ujfpApi + Q + 2A: + 36^ + 2^ + ^k^ -{e + 2k + ek^) u? + ^-u^^p,^^fl^ 
-\ (e -uj^ + 2^)pop; -2kh (uj) 6 {uj^ - e)p[ - e)podkp[ + 2upodkp[ ; (45) 



k 
and: 

/a (w, uj'] k) =Y^ 

e=± 



'--k^ - 6ek'^ - -k^ - (6ek - 6k^ + 2ek^) + (9 + k'^) + ^ekuj^ 
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\2k^ k 3 2 
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2k 
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15 2 



A;3 A; 



A;3 



+ 7 N + 7^^ P^p't + (6 - 12t^) Po^fePj 



+ ( -3 + e/ccj + —uo 



6 



po^fcp; + 3kpodlpi - ^ (^^ - Po^fePj 



+12kh (u) 5 {u^ - k^) p'l-Q \uj\ 6 (cu' - k^) p[ - 6k^ \uj\ 8^26 (lu^ - k^) p'l . 



(46) 



It remains to perform the integrals over the frequencies uj and uj' and then over the 
momentum k. Of course, these integrations are not straightforward and necessitate numerical 
work Also, the dimensionless parameter m{Nc,Nf) 



m, 



4 



IS 



implicitly present in the spectral densities pi^t and so, each case has to be treated separately. 
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Recall that this direct calculation is performed in the sole context of the Braaten-Pisarski 
HTL-summed next-to-leading order perturbation. As we emphasized in the introductory 
remarks, there is the problem of occurrence of infrared divergences to be aware of. Hence, 
extra work is needed in order to extract these from the finite contributions. One interesting 
point to wander about is what sort of divergences we will obtain. Indeed, in the direct 
calculation of 7/(0), the damping rate for longitudinal gluons with zero momentum, the 
divergent term behaves like 1/rf |^|. In the second coefficient in of 'jt{p), the damping 
rate for transverse gluons, l/rj'^ does also appear together with Inr] The question is 

then: what sort of divergences will we get for 7±(p)? If different from l/rj'^ and In 77, are we 
able to understand why? 

We stress once again that the occurrence of these divergences may simply be due to the 
early expansion we make of the HTL-summed next-to-leading order self-energies in powers 
of the external momentum. We have argued otherwise in section three, but a really more 
convincing argument would be to carry the very same calculation in a way that avoids such 
an early expansion. We have indicated that, in the perturbative context, this could be 
technically very difficult. 

In any case, there is by now convincing evidence in the literature that next-to-leading 
order quantities are magnetic-sensitive. We tend to be of the viewpoint that the occurrence 
of infrared divergences in HTL-summed next-to-leading order self-energies is probably a 
manifestation of this magnetic sensitivity, and that a more complete next-to-leading order 
calculation should remove them. We mean a calculation that takes into account magnetic 
effects not present at the electric scale, and hence not incorporated in the hard thermal 
loops. Of course, we may need to understand first such effects. It is then interesting to ask 
whether the infrared divergences one obtains, in the damping rates and in possibly similar 
quantities, can be of any help. It may also be that these effects cannot even be grasped 
perturbatively. 

Also, we have argued in section three that regularizing the infrared region with a simple 
shift in the static transverse gluonic propagator by a momentum-independent magnetic mass 
may not be the best way to shield from magnetic sensitivity. This is important in view of 
the fact that much of the results of high temperature QCD rely on such a regularization. 
In our work, we exclude from the outset the magnetic region by introducing an infrared 
cutoff of the order of the magnetic scale. One drawback is that our calculation is valid 
only for very soft external momenta and the results we obtain cannot be carried to larger 
values. That our results differ analytically from other estimations, particularly those relying 
on the magnetic-mass shielding, comes mainly from the regularization procedure and its 
implications. Also, the difference in the results accentuates the sensitivity of the HTL-based 
perturbation to the magnetic scale and is added evidence that possibly interesting physics 
is happening between the soft and very soft regions. 
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